t=-1:0.00001:1;  
A0=10;  %载波信号振幅
f=6000;  %载波信号频率
w0=f*pi;
Uc=A0*cos(w0*t);     %载波信号
figure(1);
subplot(3,2,1);
plot(t,Uc);
title('载频信号波形');
axis([0,0.01,-15,15]);
subplot(3,2,2);
Y1=fft(Uc);             %对载波信号进行傅里叶变换 
plot(abs(Y1));title('载波信号频谱');
axis([5800,6200,0,1000000]); 

t=-1:0.00001:1;  
A1=5;   %调制信号振幅
f=6000;  %载波信号频率
w0=f*pi;
mes=A1*cos(0.001*w0*t);   %调制信号
subplot(3,2,3);
plot(t,mes);
xlabel('t'),title('调制信号'); 
subplot(3,2,4);
Y2=fft(mes);             % 对调制信号进行傅里叶变换 
plot(abs(Y2));
title('调制信号频谱');
axis([198000,202000,0,1000000]);

t=-1:0.00001:1;  
A0=10;  %载波信号振幅
A1=5;   %调制信号振幅
A2=3;   %已调信号振幅
f=3000;  %载波信号频率
w0=2*f*pi;
m=0.15;    %调制度
mes=A1*cos(0.001*w0*t);   %消调制信号
Uam=A2*(1+m*mes).*cos((w0).*t); %AM 已调信号  
subplot(3,2,5);
plot(t,Uam);
grid on;
title('AM调制信号波形');
subplot(3,2,6);
Y3=fft(Uam);              %  对AM已调信号进行傅里叶变换  
plot(abs(Y3)),grid;
title('AM调制信号频谱');
axis([5950,6050,0,500000]); 
